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1 .0  INTRODUCTION 


Our  knowledge  of  unsteady  boundary  layers  has  been  enhanced  In  recent  years  by 
a  combination  of  experimental  and  computational  Investigations.  The  former 
have  been  concerned  mainly  with  the  problem  of  airfoils  which  oscillate  In 
subsonic  flows  but  Include  a  small  number  of  transonic  flow  Investigations. 

The  latter  have  Involved  the  solution  of  boundary-layer  equations  for  flows 
over  cylinders  and  airfoils  and  have  been  concerned  largely  with  the  nature  of 
the  singularity  which  occurs  near  separation.  Solutions  of  the  Navier-Stokes 
equations  have  also  been  obtained  to  provide  an  overall  understanding  of  the 
flow  patterns  at  conditions  near  or  at  dynamic  stall.  Useful  reviews  of  prev¬ 
ious  work  have  been  provided  by  Telionls  (1979),  McCroskey  (1982),  Williams 
(1977)  and  Cebecl  (1982). 

The  experimental  Investigations  have  been  carried  out  mainly  to  Improve 
understanding  of  flow  around  helicopter  rotor  blades  and  include  the  low-speed 
measurements  of  McCroskey  et  al.  (1982),  Carr  et  al.  (1977),  Carr  and  McAlister 
(1983),  Young  (1982),  Geissler  (1983)  and  Coustelx  et  al.  (1981)  and  transonic 
flow  measurements  of  Tijdeman  (1977),  Davis  and  Malcolm  (1979).  The  range  of 
measurements  of  McCroskey  et  al.  and  Carr  et  al.,  in  particular,  is  extensive 
and  Includes  more  than  fifty  combinations  of  Mach  number  and  parameters  of  the 
unsteady  motion  for  each  of  eight  airfoil  sections.  As  a  consequence,  four 
flow  regimes  have  been  Identified  and  correspond  to  no  stall,  stall  onset, 
light  stall  and  deep  stall.  It  appears  that  the  breakdown  of  the  unsteady 
boundary  layer  leads  to  a  large  vortex  which  Is  formed  near  the  surface  at 
large  angles  of  attack  and  causes  stall  to  occur  shortly  thereafter. 

Computational  Investigations  of  oscillating  airfoil  flows  have  been  of  two 
main  types.  In  the  first,  solutions  of  the  Navier-Stokes  equations  have  been 
obtained,  for  example  by  Mehta  (1977)  for  Incompressible  laminar  flows  and  by 
Shamroth  (1981)  for  compressible  turbulent  flows.  Further  studies  of  this 
type  are  clearly  necessary  and  results  so  far  show  qualitative  features  of  the 
flow  field  and  lift  curves  which  are  again  in  qualitative  agreement  with 
experiment.  The  second  approach  has  Involved  the  solution  of  the  boundary- 
layer  equations.  The  Investigations  of  Cebecl  and  Carr  (1981,  1983)  have 
reported  results  for  an  external  velocity  distribution  typical  of  those  found 
near  the  leading  edge  of  thin  airfoils.  The  existence  of  a  singularity  In  the 


solutions  Is  evident  In  the  vicinity  of  the  leading  edge  at  the  higher  angles 
of  attack  and  has  led  to  the  subsequent  Investigations  of  Its  nature,  see  for 
example,  Cebecl,  et  al.  (1983). 

The  need  for  more  fundamental  Investigations  of  the  use  of  boundary- layer 
equations  to  represent  unsteady  flows  Is  clear  from  the  airfoil  Investigations 
of  the  previous  paragraph.  In  this  connection,  the  oscillating  flat  plate 
Investigations  of  Coustelx  et  al.  (1981)  and  the  contributions  of  Tellonls 
(1974),  van  Oommelen  and  Shen  (1982a, b).  Smith  (1982)  and  Cebecl  (1979,  1982) 
for  flow  over  a  cylinder  Impulsively  started  from  rest  are  particularly 
relevant.  In  addition,  Williams  (1982)  and  Williams  and  Stewartson  (1983) 
have  made  Important  contributions  to  our  understanding  of  the  nature  of  the 
singularity  and  Its  consequences.  Perhaps  the  most  Important  contribution  has 
been  that  of  van  Dommelen  and  Shen  (1982a)  who  solved  the  boundary- layer 
equations  In  Lagranglan  form  and  revealed  the  existence  of  the  singularity. 

This  result  suggests  the  need  for  Interaction  between  the  viscous  and  Invlscld 
flow  equations. 

Though  a  number  of  studies  have  been  conducted  to  Improve  our  understanding  of 
unsteady  two-dimensional  boundary  layers,  very  little  work  has  been  done  for 
unsteady  three-dimensional  flows.  Experimental  Information  Is  lacking,  but 
the  similarity  between  the  two-  and  three-dimensional  equations  suggests  that 
the  same  phenomena  may  occur.  The  three-dimensional  boundary  layer  on  a  body 
of  revolution  at  angle  of  attack  Is  obviously  more  complicated  than  for  air¬ 
foils  or  cylinders  but  Is  similar  on  the  line  of  symmetry,  provided  that  cross- 
flow  gradients  are  taken  Into  account.  In  the  case  of  steady  flow  past  slender 
spheroids.  It  Is  known  [see,  for  example,  Cebecl  et  al.  1980]  that  separation 
does  not  occur  on  the  leeward  side  until  well  past  the  maximum  thickness  If 

the  angle  of  attack  a  Is  less  than  a  critical  value  o  (*41°),  but  for  a  >  a 

c  c 

there  Is  a  dramatic  change  with  separation  occurring  very  close  to  the  nose. 

If  the  external  velocity  Is  prescribed,  this  abruptly  terminates  the 
Integration.  For  thin  airfoils  there  Is  a  parallel  situation,  as  shown  by 
Cebecl  et  al  (1981a),  but  the  critical  angle  Is  now  of  the  order  of  magnitude 
of  the  airfoil  thickness. 

A  study  of  the  unsteady  boundary  layer  on  the  line  of  symnetry  to  determine 
the  relation  between  unsteady  separation  and  singularities  In  the  solution  Is 
reported  here.  For  the  first  step  we  concentrate  on  the  separation 
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problem  and  do  not  consider  the  effect  of  moving  stagnation  points  which  Is 
Important  but  adds  considerably  to  the  computational  difficulties.  The  par¬ 
ticular  problem  we  study  Is  the  development  of  a  boundary  layer  on  a  thin 
spheroid  In  uniform  motioh  at  constant  angle  of  attack  after  an  Impulsive 
start.  Angles  of  attack  ranging  from  30°  to  50°  are  chosen  for  this  purpose. 

Of  these  a  *  45°  corresponds  to  nose  separation  which  Is  marginal  while  a  *  50° 
corresponds  to  a  strong  steady-state  singularity.  The  problem  has  been  formu¬ 
lated  for  the  general  case  of  spheroids  at  angles  of  attack  and  specialized  to 
thin  bodies,  and  specifically  those  for  which  the  thickness  ratio  t  Is  van¬ 
ishingly  small. 


2.0  APPROACH 


With  respect  to  a  cylindrical-polar  set  of  coordinates  (x,r,e)  and  origin  0, 
let  the  equation  of  the  spheroid  be 


where  x  is  a  positive  number  less  than  unity,  which  may  be  regarded  as  a 
measure  of  the  thickness  ratio  of  the  prolate  spheroid.  Here  we  have  nondimen* 
sionalized  the  coordinates  and  time  t  taking  the  semi-major  axis  and  the 
velocity  at  infinity  to  be  unity.  We  assume  that  at  an  infinite  distance  from 
the  spheroid  the  velocity  of  the  fluid  is  in  a  direction  lying  in  the  merid¬ 
ional  plane  and  at  an  angle  a  to  Its  major  axis.  Further  we  assume  that 
outside  the  boundary  layer  the  Invlscld  flow  Irrotational  and  attached  and 
neglect  the  circulation  around  the  spheroid.  The  achievement  of  our  goal  of 
determining  the  development  of  the  boundary  layer,  given  the  main-stream 
velocity,  would  be  an  Important  step  forward  in  the  task  of  predicting  the 
circulation  around  bodies  of  revolution  at  Incidence. 

With  these  assumptions  the  velocity  of  slip  on  the  spheroid,  according  to 
invlscld  theory,  has  components 

ue  =  V0(x)  coso  cosg  -  Vgo(x)  sino  sine  cose 

we  *  Vgo(x)  sino  sine 


where 


(2) 


Figure  1.  Notation  for  prolate  spheroid. 

and  a  is  a  specified  function  of  time,  being  constant  in  the  present  study. 

The  configuration  in  the  meridional  plane  Is  illustrated  in  figure  1  with  S  the 
stagnation  point  of  the  inviscid  flow  defined  by  e  *  0,  x  *  -{1  -  x^  tan^a}. 


tm 
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Turning  to  the  boundary  layer  we  define  (u,v/v,w)  to  be  the  velocity  components 
respectively  along  the  meridional  lines  0  *  const,  along  the  normals  to  the 
spheroid  and  in  the  azimuthal  direction  where  v  is  the  kinematic  viscosity  of 
®  the  fluid  and  is  assumed  small.  Then 


lx  (h2u)  +  55  (hl"1  *  5y  (l,lh2vl  ■  0 

*  9 

3u  ,  u  3U  .  W  3u  .  „  3U  .  ylu,  1  3p  .  3  U  , 

7T  +  RJ  3x  +  TiJ  55"  +  v  ?y  +  "^2  +  J  1 

*  9 

3U  U  3W  W  3W  3W  m  1  3£  .  3  W  , 

3T*K75x*K7ae  +  vW'u*Z  ^  1 

where  y/v  measures  distance  along  the  outward-drawn  normal  from  the  body. 


Is  the  geodesic  curvature  of  the  surface  lines  x  =  const,  and  p  Is  the  reduced 
pressure.  The  appropriate  boundary  conditions  are 


u  =  v  =  w  =  0  at 

u  ue,  w  Wg  as 


y  *  0 
y  -*•  ® 


(10a) 


The  specification  of  the  problem  Is  completed  by  assigning  initial  conditions 
on  u,w  with  respect  to  space  and  time.  The  spatial  conditions  are  that  on  the 
line  of  symmetry  w  «  0  and  u  vanishes  for  all  y  at  that  value  of  x  for  which 
ue  *  0.  This  is  physically  and  mathematically  sensible  for  impulsive  prob¬ 
lems  in  which  the  stagnation  point  is  fixed  but  is  not  necessarily  correct 
when  the  angle  of  attack  is  varying  with  time,  see  Cebeci  and  Carr  (1981)  for 
example.  The  temporal  conditions  are  that  at  t  -  0  both  u  and  w  vanish  on  the 
body  y  *  0  on  the  line  of  symmetry  but  are  equal  to  their  external  values 
elsewhere.  Thus 


u  =  w  *  0  at  •  y  ■  0,  t  »  0  for  all  x 


while 


(10b) 


u  *  ug,  w  *  wfi  for  y  >  0,  t  »  0  and  all  x. 

In  the  early  studies  of  the  properties  of  the  boundary  layers  on  prolate 
spheroids  conducted  by  Wang  (1970)  and  Hirsh  and  Cebeci  (1977),  some  diffi¬ 
culty  was  experienced  in  continuing  the  solution  past  the  nose  at  x  =  -1,  pri¬ 
marily  because  of  the  singularities  In  the  properties  of  hj,  h2  and  K2  there. 

A  common  procedure  was  first  to  perform  the  Integration  along  the  line  of  sym¬ 
metry  from  the  stagnation  point  to  as  near  the  nose  as  possible  and  then  to 
jump  across  the  nose  to  the  same  value  of  x  on  the  leeslde  (9  ■  m)  assuming 
that  the  flow  properties  are  essentially  unchanged.  After  that  Integration  on 
the  leeslde  may  be  continued  as  far  as  separation.  Afterwards  the  procedure 
may  be  extended  to  more  general  points  In  the  neighborhood  of  the  nose.  This 
technique  Is  effective  at  moderate  values  of  t  [t  ■  1  corresponds  to  a  sphere] 
but  leads  to  difficulties  as  t  ♦  0  especially  at  high  angles  of  attack.  Cebeci 
et  al.  (1980)  demonstrated  that  the  singularity  may  be  removed  by  a  suitable 


transformation  of  the  surface  coordinate  system  enabling  a  smooth  passage  to 
be  made  around  the  nose.  The  transformation  Is  equally  effective  for  unsteady 
boundary  layers  and  Is  now  explained  for  the  limiting  case  of  a  paraboloid 
which  corresponds  to  x  ■  0. 

The  first  step  Is  to  define  new  surface  coordinates  by 

x2-l 


2  t(1  -  x") 

(11 

X  =  S  cose,  Z  =  S  sine,  Y  =  y/x,  T  =  t/x2 
with  5  =  0  at  x  =  -1 ,  and  new  velocity  components  by 

u  *  u  cose  +  w  sine,  w  *  w  cose  -  u  sine,  v  =  v/x  (12 

The  purpose  of  this  transformation  is  to  convert  the  polar  form  of  the  equa 
tlons  (5)  -  (7)  into  a  quasl-rectangular-Cartesian  form  which  is  free  of 
singularities.  The  governing  equations  reduce  to 

N(|y  +  fy)  +  |y-L(UX  +  WZ)-0  (13 

|«.N(U|U  +  u||)tLW(Wx.uZ)t,|«.^+A  ,H 

!¥  +  N(u|y  +  -  UHWX  -  UZ)  +  V  IS  ■  (15 


where  ^  ,  &  are  pressure-gradient  parameters  defined  by  the  requirement 
that  u  1  ue?  w  +  we  as  Y  >  ®  and 

N  *  ’  L  s  5"  ^ +  K2)  (16 

For  the  paraboloid  we  let  x  -*■  0  after  defining 

(1  -  x2)^2 

P  ■  -  H7 


exp  C(1  +  p2)1/2  -  1] 


(18 


N 


,  _  (1  +  P2)V2  -  1 

L  -  '2J/2 


psn  +  p‘r 


(19) 


and  so  Eqs.  (13)  -  (15)  are  independent  of  x.  We  write  (Ug,  Wg)  for  the 
limits  of  (ug  cose  -  wg  sine,  ug  sine  +  wg  cos  e)  as  t  ♦  0  and  then 

U  »  PX  cos%  T/-  -  2(1  -  5^)  sinct  (20) 

e  S(1  +  p2)1/Z  * 


We 


pZ  C0Sa  .  2XZLp  ... 
r  ?  I/?  +  — sma 

S(1  +  p 2)1/2  S 


(21) 


so  that 


11  /  it 


3U. 


3U 


(22) 


3W  3W 

$2  s  N*ue  TT  +  we  ~  LMwex  "  uez^  ^23* 

These  equations  are  free  of  all  geometric  singularities  and  In  particular  at 
the  nose  p  =  0.  There  Is  therefore  now  no  special  problems  about  integrat¬ 
ing  the  equations  through  the  nose  although  it  should  be  noted  that  the  equa¬ 
tions  are  only  appropriate  at  distances  from  it  corresponding  to  p  *  0(1). 

This  is  the  natural  length  scale  for  the  paraboloid  but  in  terms  of  general 

2 

thin  axi symmetric  bodies  this  distance  is  0(t  )  from  the  nose. 


2.1  Line-of -Symmetry  Equations 

These  may  be  deduced  from  the  general  form  Eq.  (13)  et  seq.  given  In  the 
previous  section  using  similar  arguments  to  those  in  Cebeci  et  al.  (1981).  We 
shall  not  repeat  them  in  detail  but  merely  state  the  transformations  and  the 


final  form.  First  we  write 

U  -  U0(X,Y,T)  +  0(Z2), 

V  *  V0(X,Y,T)  +  0(Z2) 

(24) 

W  -  Z  exp[l  -  /I  +  p^]Wj(X,Y,T)  +  0(Z3) 

(25) 

and  allow  for  negative  values  of  X  by  permitting  p  to  become  negative.  When 
p  <  0  the  sign  of  S  In  Eq.  (18)  must  be  changed  and  generally  X  ■  S  sgnp  In 
the  limit  Z  ♦  0.  On  the  leeslde  line  of  symetry  p  <  0. 


8 


We  now  substitute  Eqs.  (24),  (25)  into  (13)  et  seq  and  take  the  limit 
Z  0.  The  external  velocity  components  reduce  to 


cosa  -  2  sina 


(1  +  p£) 


u  _  [(1  +  p2)^2  +  l]cosa  +  2p  sina 

’•  (TT7?77 


To  put  the  equations  into  the  most  convenient  form  for  numerical  Integra-  tion 
of  the  impulsive  problem,  with  primes  denoting  differentiation  with  respect  to 
n  we  write 


(1  ♦  p  ‘)"V 


»  T  =  T"1/2(l  +  p2)1/4V  -  ueT 

ti  r  rt 


UosV'  ’ 

e 


0  2(1  +  p2) 

Wt  *  g* 


where  f  and  g  are  functions  of  p.n.t.  Then 


\  = '  V  !£  • f {alV  +  “3jT  ]  '  a29 


where 


'  2(1  +  p2)  (p2  +  1),/Z  *  1  ’  **  1  +  (lYp2) 

and  the  appropriate  boundary  conditions  for  f  and  g  are: 


f  •  f  -  g 

» g' 

=  0 

at 

J3 

II 

O 

• 

f'  -  i, 

g* 

*  Wle 

as 

n  -*•  00 

The  momentum  equations  reduce  to 

f-  ♦  P,ff  *  P6gf>  *  P2[l  -  (f)2]  p  P„nf  -  P7  |f  p  P,0(f  ff  -  f  ff 

(32) 

g-'  +  P^fg"  +  P6gg-  +  P3(f')2  +  P4(g')2  +  P5f'g'  +  PT2  +  P13ng" 
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These  two  equations  (32)  and  (33)  are  In  a  convenient  form  for  numerical 
Integration  since  at  T  =  0  f,g  and  their  derivatives  are  given  by 

2  2 

f  =  nerf(Xn)  -Ul  -  5X  n  ),  g  -  W,  f 

”  X/tr  le 

f  =  erf(xn),  g'  *  Wlef  (35) 

,2  2 

f «  =  e’x  n  ,  g-  .  W,  f- 

/*  le 

where 


X  .  J  (1  ♦  p2)1" 

Further,  the  Integration  In  the  direction  p  Increasing  may  be  assumed  to  start 
at  the  stagnation  point  p  ■  2  tana,  where  f,g  are  Independent  of  p  and  sat¬ 
isfy  partial-differential  equations  In  n*T  only,  and  proceeds  In  a  straight¬ 
forward  way  since  UQ  >  0  when  p  >  2  tana.  A  similar  remark  applies  to  the 
Integration  In  the  direction  of  p  decreasing  since  U  <  0  when  p  <  2  tana. 

''a 
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3.0  SOLUTION  PROCEDURE 


The  numerical  solution  of  equations  (32)  and  (33)  was  obtained  using  the  Box 
method,  a  two-point  finite-difference  method  developed  by  H.  B.  Keller  and 
extensively  used  for  two-dimensional  flows,  both  steady  and  unsteady,  and  for 
three-dimensional  steady  flows.  Here  we  use  two  versions  of  the  Box  method, 
the  regular  box  In  regions  where  f 1  _>  0  across  the  layer  and  the  zig-zag  box 
in  regions  where  f'  <  0  at  any  value  of  r\.  Specifically  in  advancing  the 
solution  to  a  new  value  of  T  or  p,  this  property  of  f*  determined  the  choice 
of  the  box  scheme.  The  details  of  the  box  scheme  are  described  in  Bradshaw  et 
al.  (1981). 

Calculations  of  the  unsteady  boundary  layers  were  carried  out  for  angles  of 
attack  a  of  30°,  40°,  45°,  50°.  Provided  flow  reversal  did  not  occur,  the 
solution  of  the  equations  by  the  regular  box  was  appropriate,  no  numerical 
difficulties  were  encountered  and  the  solution  was  smooth.  The  step-sizes 
chosen  for  all  computations  were  Ap  »  0.2,  AT  =  0.2;  in  selected  cases  compar¬ 
isons  were  made  with  studies  using  step-sizes  of  Ap  =  0.1  and  AT  =  0.1  or 
smaller  and  the  differences  were  negligible.  However,  once  flow  reversal 
occurred,  the  regular  box  became  prone  to  instabilities  and  was  replaced  in 
this  domain  by  the  zig-zag  box  over  all  n  and  the  step-sizes  in  p  and  T  were 
reduced.  Various  calculations  were  carried  out  with  a  nonuniform  mesh  in  p,  T 
in  the  neighborhood  of  separation.  The  smallest  value  of  Ap  was  0.0005  near 
the  separation  point;  the  step  size  in  T  was  progressively  reduced  as  the 
calculations  proceeded,  the  smallest  value  being  0.01.  Elsewhere  Ap  was 
fixed  at  0.2.  Comparisons  with  solutions  using  larger  step-sizes  gave 
confidence  that  the  results  were  reliable  outside  the  separated  region  for 
T  <  7,  but  that  inside  the  separated  region  smaller  step-sizes  In  p  and  T 
were  necessary  to  avoid  small  oscillations  which  developed  for  T  >  5.5  at 
a  »  50°. 


4.0  RESULTS 


The  principal  properties  of  the  unsteady  boundary  layers  are  the  skin- 
friction  components  defined  by 
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Uq(o,P,T)  =  T‘1/2Ue  f"(o,p,T)  =  (1  +  p2)1/4  (^)w 

0 

(36a) 

W ( o , p , T )  =  T~1/2g"(o,p,T)  =  (1  +  p2)1/4  (-^-)w 

(36b) 

and  the 

displacement  thicknesses  defined  by 

A  -|(p,T) 
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These  are  displayed  in  Fig.  2  for  a  =  30°  when  separation  does  not  occur. 

I 

In  Fig.  3  we  display  the  properties  of  UQ  for  a  *  40°,  a  solution  which 
is  on  the  verge  of  separation,  and  in  Fig.  4  we  show  all  the  functions  [Eqs. 

(36)  and  (37)]  for  a  *  50°  when  separation  clearly  occurs.  Where  appropriate, 
the  corresponding  steady-state  results  of  Cebeci  et  al.  (1980)  are  included  in 
these  graphs  for  comparison  purposes. 

The  graphs  show  that  for  p  >  0,  l.e.  on  the  windward  line  of  symmetry,  the 
steady-state  solution  is  rapidly  approached  as  T  increases  and  is  essen¬ 
tially  established  for  T  >  1. 

For  p  <  0,  i.e.  on  the  leeward  line  of  symmetry,  the  flow  properties  In  the 
steady  state  have  two  interesting  features.  First  for  any  value  of  a,  the 
cross-flow  displacement  thickness  a2  rapidly  decreases  with  p  and  soon 
becomes  negative,  even  If  separation  has  not  occurred,  and  In  fact  Is  exponen¬ 
tially  large  and  negative  when  p  Is  large  and  negative.  The  unsteady  data  for 
a  ■  30°,  Fig.  2d,  confirms  this  trend  but  the  limits  T  finite  and  p  -►  -®  for 
the  unsteady  flow  and  p  ®  for  the  steady  flow  are  not  Interchangeable.  The 
steady-state  solution  for  a2  becomes  exponentially  large  and  negative  for 
p  >  -2,  and  the  approach  of  the  unsteady  solution  to  the  steady  state  has  the 
unusual  feature  In  that  for  some  T,  say  T  >  2,  the  steady  and  unsteady  solutions 


— 


agree  quite  well  up  to  some  negative  p  after  which  the  unsteady  solutions 
flatten  out  as  they  approach  p  -*•  -  «.  As  T  Increases,  the  divergence  of  steady 
and  unsteady  solutions  occur  at  progressively  more  negative  Ag,  and  as  T  ® 
the  steady  and  unsteady  solutions  agree  up  to  a2  *  -00. 


Second  the  steady-state  meridional  component  UQ  of  the  skin  friction  does  not 
tend  monotonlcally  to  a  limit  as  p>  -®,  but  even  at  small  values  of  a  has  a 
trough  followed  by  a  peak  due  to  the  overshoot  In  the  external  velocity  (26). 


Further,  at  a  >  a  *  41°,  U  actually  vanishes  at  p  *  p  (a),  when  p 

v  O  C  v 

Is  a  point  near  p  *  -1  depending  on  the  particular  choice  of  a.  The  solution  is 
singular  at  this  point  and  the  calculation  must  terminate. 


The  unsteady  solution  approaches  the  steady  state  as  T  ®  for  a  <  ac  (Figs.  2 
and  3)  but  more  slowly  on  the  windward  side.  Some  overshoot  may  occur  but  it 
Is  small  and  may  be  an  Indication  of  a  fall-off  In  accuracy  of  either  the 
steady-state  or  the  unsteady  solutions.  For  a  >  a  (Fig.  4),  the  same 


remarks  apply  as  long  as  p  >  p  (a)  and  flow  reversal  has  not  occurred.  For 

c 


smaller  values  of  p,  reversed  flow  occurs  and  there  Is  no  corresponding  steady- 
state  solution.  If  a  »  45°,  separation  first  occurs  at  T  *  5.64  when  p  *  -1.2; 
It  remains  confined  within  the  range  -1.0  >  p  >  -1.4  and  Is  weak  until  T  «  6 

I 

when  the  computation  was  terminated.  It  might  be  argued  that  W-j  and  a2  are 
approaching  limit  states  as  T  »  but  A^  Is  definitely  showing  signs  of  a  pro¬ 
nounced  negative  minimum  near  p  »  -1.2. 


Separation  is  marginal  for  a  »  45°  and  so  a  more  extensive  computation  was 
carried  out  for  a  *50°.  In  this  case  separation  set  in  when  T  *  3.71  at 
p  *  -1.083  and  gradually  spread  out  as  T  Increased  to  extend  around  the  range 
-0.908  £  -1.396  at  T  «  6  when  the  calculations  became  somewhat  dubious 

due  to  oscillations  and  Instabilities,  a  behavior  previously  observed  In 
relation  to  the  circular  cylinder  started  Inpulslvely  from  rest,  Cebecl 

I  . 

(1982).  The  behavior  of  the  separated  solutions  of  UQ  for  p  <  -1  and  4  < 

T  £  6  are  shown  In  figures  4a  and  5b,  Indicating  that  reattachment  occurs  In  this 
time  Interval.  It  is  Inferred  that  after  slightly  larger  time  than  T  ■  6, 
reattachment  does  not  occur  and  separation  approaches  that  of  steady  state 
(p  *  0.91),  as  corroborated  by  figure  5a  where  Aj  develops  a  rapid  growth  rate. 
The  curves  of  a2,  not  shown  for  p  <  -1  In  figure  4d,  do  not  develop  any  unusual 


a 
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features,  they  become  monotonies! 1y  more  negative  as  p  decreases,  and  In  addl 
tlon  the  curves  for  T  >  4  are  practically  conlcldent  with  the  T  »  4  curve  to 
at  least  p  «  -2.  Over  the  range  of  T  considered  (T  <  6)  the  variation  of 
and  appears  to  remain  smooth  but  UQ  and  Aj  develop  sharp  minima  just 
downstream  of  separation. 


5.0  THE  SINGULARITY? 


We  now  Investigate  whether  the  results  obtained  for  a  ■  50°  allow  us  to 
conclude  that  when  the  computation  terminated  the  flow  properties  are  consist¬ 
ent  with  a  singularity  centered  just  downstream  of  separation  In  the  neighbor¬ 
hood  of  T  *  6. 


Two  possible  structures  for  a  terminating  singularity  to  the  solution  have 
been  proposed  of  which  one  Is  due  to  Smith  (1982).  In  his  model,  the  viscous 
terms  are  significant  at  breakdown,  which  Is  a  generalization  of  the  Goldstein 
singularity,  and  both  skin-friction  and  displacement  thickness  are  Infinite. 

No  example  of  such  a  structure  has  yet  been  found  and  although  we  cannot  rule 
out  its  relevance  In  the  present  flow,  our  Inclination  Is  to  prefer  the  second 
possibility,  due  to  van  Donmelen  and  Shen  (1982a),  In  which  the  skin  friction 
remains  finite  while  the  displacement  thickness  becomes  Infinite.  Strong 
evidence  In  favor  of  the  unsteady  boundary  layer  on  a  circular  cylinder 
behaving  In  the  same  way  has  been  provided  by  van  Dommelen  and  Shen,  Cowley 
(1983)  and  Ingham  (1984)  and  on  an  oscillating  airfoil  by  Cebecl  et  al.  (1983) 


The  essence  of  the  singularity  Is  Invlscld  and  Independent  of  the  pressure 
gradient.  Thus,  In  the  case  of  a  two-dimensional  boundary  layer,  with 
equations 
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let  us  suppose  the  singularity  Is  centered  at  t  »  t  ,  x  *  xQ.  Then  In  Its 
simplest  form  we  write 

u  -  -u0  +  (^  -  t)V2  u(x,y)  (39) 

when  0  <  t0  -  t  «  1,  |x  -  x0|  «  1, 

where 


x  «  [(x  -  x0)  -  u0(t0  -  t)]/(t0  -  t)3/2,  y  »  y(t0  -  t)V4 


and  the  determination  of 'b  reduces  to  quadratures.  This  structure  must  be 
matched  to  other  structures  near  y  ■  0  and  as  y  ♦  ®  but  the  principal  gross 
features  can  be  Inferred  from  (39).  The  displacement  thickness  has  a  peak  at 

'x*  -7/4 

x  *  0  which  becomes  Infinite  like  (tQ  -  t)  '  as  t  tQ.  On  the  other  hand, 
the  skin  friction  remains  finite  as  t  -*•  t0  for  all  x.  The  velocity  profile 

becomes  very  flat  as  t  -*•  t  provided  x  1,  with  u  ■  -u„  over  a  length 

-1/4  0  0 

(tQ  -  t)  '  normal  to  the  surface,  the  adjustment  to  the  boundary  conditions  both 

on  the  surface  and  at  Infinity  taking  place  over  effectively  finite  ranges  of 

values  of  y.  Further  details  of  the  structure  of  the  singularity,  on  the  basis 

of  a  Lagrangian  representation  may  be  found  In  van  Dommelen's  thesis  (1981) 

and  in  Elliot,  et  al.  (1983)  on  the  basis  of  an  Eulerian  representation  as 

exemplied  by  (38). 

This  theory  may  be  adapted  to  our  problem  at  least  in  principle.  We 
replace  u,x,y  In  (39)  by  f',p,n  and  neglect  the  left-hand  side  of  (32)  in 
the  neighborhood  of  p  *  p0,  t  =  t0.  Then  the  equation  for  is  the  same 
as  that  for  u  and  we  may  infer  that  f*  Is  as  likely  to  develop  the  van 

Dommelen  and  Shen  singularity  as  is  u  in  the  circular  cylinder  problem.  The 

displacement  thickness  Aj  can  be  expected  to  develop  an  Increasingly  sharp 
peak  as  t  is  increased  whereas  the  corresponding  component  of  skin  friction 
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U  remains  finite  and  smooth, 
o 

The  evidence  in  Fig.  4  supports  these  conclusions,  but  a  more  severe  test 
Is  to  plot  the  streamwlse  displacement  velocity,  3Aj/3p,  as  a  function  of  p  for 
various  T  and  this  Is  done  in  Fig.  5a  after  smoothing  the  data  for  p  <  -0.95. 

The  results  of  Cowley  (1983)  for  the  solution  of  (38)  with  ug  *  sinx  show  that 
the  singularity  occurs  at  T  *  3.0.  The  similarities  in  form  are  striking;  the 
principal  difference  being  a  slight  drift  In  3Aj/3p  away  from  separation  as  T 

Increases.  The  more  detailed  plots  of  in  Fig.  5b  confirm  the  smooth  behavior 

of  this  function. 

If  now  we  examine  (33)  we  see  that  the  same  form  for  g'  may  be  assumed 

I  I 

provided  gQ  Is  taken  equal  to  fQ  and  again  the  left-hand  side  of  (33)  may  be 
neglected.  This  means  that  a2  should  also  develop  a  needle-like  singularity  and 
Fig.  4d  shows  no  evidence  of  this.  However,  plots  of  3A2/3p  as  functions  of  p 
for  various  T,  displayed  In  Fig.  5c,  show  clear  signs  of  the  development  of  a 
strong  minimum,  similar  to  that  of  Cowley  (1983),  which  does  drift  towards 


separation.  The  alternative  possibility  that  Lz  1s  finite  seems  to  require  t! 

as  n  increases  g*  has  reached  the  value  Wj  before  the  inviscid  zone  n  -  1  is 
reached.  Were  this  the  case,  the  two  sides  of  (33)  would  not  balance  over  th 

I 

since  then  f'  “  f.  f  1. 


6.0  DISCUSSION 


In  this  report  we  have  examined  the  evolution  of  the  unsteady  boundary  layer  on 
the  line  of  symmetry  of  a  paraboloid  which  Is  set.  Impulsively,  Into  motion  at 
T  «  0  with  uniform  velocity  and  at  an  angle  a  of  attack.  We  have  shown  that 
If  a  <  a  (=41°)  and  the  steady-state  boundary  layer  exists  everywhere  on  this 

C 

line,  the  unsteady  solution  approaches  It  reasonably  quickly  and  without  any 
significant  special  features.  The  same  Is  true  for  a  >  ac  on  the  windward 
side.  On  the  leeward  side,  the  steady-state  boundary  layer  separates  and  the 
solution  must  terminate  there  because  the  external  pressure  gradient  Is  fixed. 
The  unsteady  boundary  layer  Is  Initially  unseparated  but  develops  a  region  of 
reversed  flow  after  a  finite  time.  A  short  time  later  the  first  displacement 
thickness  A-j  develops  a  pronounced  peak  and  we  advance  arguments  for 
believing  that  this  Is  associated  with  an  incipient  singularity  which  termin¬ 
ates  the  calculation. 

The  results  obtained  here  are  also  relevant  to  other  bodies,  most  obviously 
prolate  spheroids  but  others  as  well.  Our  earlier  studies  demonstrated  that 
In  the  neighborhood  of  the  nose,  the  governing  equations  for  a  spheroid  can  be 
transformed  Into  a  form  In  which  the  paraboloid  equations  appear  as  a  natural 
limit  (t  0,  see  Eq.  (17))  In  which  the  various  parameters  of  the  flow 
remain  finite.  Hence  we  may  confidently  expect  that  any  special  feature,  such 
as  separation,  which  appears  In  the  one  will  also  appear  In  the  other.  Wang 
and  Fan  (1982)  have  studied  the  unsteady  boundary  on  a  prolate  spheroid  for 
a  *  45®  for  a  value  of  t  of  0.25  and  found  that  the  unsteady  boundary 
layer  did  not  separate  for  any  t  but  that  the  steady-state  boundary  layer 
separated  at  p  «  0.39/t.  This  result  may  be  In  conflict  with  ours  and 
moreover  a  matter  of  surprise.  For  not  only  does  a  small  nonzero  x  make 
very  little  difference  to  our  form  of  the  equations,  but  the  result  Is  hard  to 
understand  In  the  general  context  of  boundary- layer  theory.  At  present  the 
discrepancy  Is  Inexplicable.  The  choice  of  coordinate  systems  for  Integration 
around  the  nose  does  seem  to  be  very  cumbersome  and  It  would  be  Interesting  to 
repeat  the  calculations  using  the  system  advocated  here.  For  the  purposes  of 
this  discussion  we  shall  set  this  calculation  at  one  side. 

We  repeat,  therefore,  that  In  our  view  the  phenomena  we  have  discussed  In 
connection  with  the  paraboloid  would  also  appear  In  a  comparable  study  of 


prolate  spheroids.  In  particular  the  singularity  would  appear  after  a  finite 
time  In  an  Integration  of  the  unsteady  equations  on  the  line  of  symmetry  In 
the  neighborhood  of  the  nose  and  just  downstream  of  separation  If  a  >  ac. 
generalization  of  this  cohcluslon  to  the  unsteady  boundary  layer  over  the  whole 
spheroid  Is  of  Interest.  In  principle  the  Integration  of  four  dimensional 
boundary- layer  equations  (x,d,y,t)  wave like  with  respect  to  (x,e,t)  using 
the  Keller-box  scheme  Is  no  more  difficult  than  those  for  steady  two- 
dimensional  flow,  provided  we  make  appropriate  modifications  to  account  for 
the  varying  direction  of  flow  across  the  boundary  layer.  It  just  takes  longer 
because  there  are  more  mesh  points  to  consider  In  the  (x,e,t)  space.  We  may 
expect  that  the  solution  of  these  equations  will  develop  a  singularity  on  the 
normals  to  a  curve  C  on  the  spheroid  after  a  finite  time,  the  precise  time 
varying  from  point  to  point  of  C,  and  it  will  be  Invlscid  In  character.  Indeed 
van  Dommelen  and  Shen  (1982b)  (see  also  Cowley  1983)  have  suggested  a  general¬ 
ization  of  their  two-dimensional  structures  which  seems  likely  to  be  appropri¬ 
ate.  From  our  experience  with  two-dimensional  unsteady  boundary  layers  we  are 
confident  that  the  Keller-box  scheme  can  carry  the  Integration  up  to  the  first 
onset  of  the  singularity.  It  remains  to  be  seen  whether  the  calculation  can 
be  extended  beyond  this  time  and  over  what  part  of  the  flowfleld  -  certainly 
the  front  portion  of  the  body  can  be  treated  but  the  precise  boundaries  are 
presumably  fixed  by  an  application  of  the  Raetz  (1957)  theory  of  Influence 
regions. 

Finally,  It  is  of  Interest  to  consider  the  possible  impact  of  this  study  on 
the  calculation  of  steady  boundary  layers  on  bodies  of  revolution.  The 
extensive  computations  In  t  on  prolate  spheroids  showed  that  on  the  windward 
part  of  the  separation  line  the  flowfleld  develops  a  Goldstein-Brown  (1965) 
singularity.  On  the  leeslde  of  the  ok  of  accessibility  (Cebecl  et  al.  1981b), 
the  calculations  for  large  a  must.  In  general,  be  terminated  by  the  external 
streamline  through  the  ok,  so  that  the  question  of  the  structures  of  the  lee¬ 
slde  separation  line  Is  Irrelevant.  At  small  values  of  a  (*6°)  this  limi¬ 
tation  does  not  apply  but  nevertheless  the  computation  breaks  down.  It  Is 
believed  that  this  breakdown  occurs  very  near  leeslde  separation  but  a  recon¬ 
ciliation  with  the  Goldstein-Brown  structure  was  not  achieved.  In  particular, 
the  blowing  velocity  Is  negative  whereas  the  theory  requires  It  to  be  large 
positive. 


A  feature  of  the  singularity  is  that  when  tQ  -  t  is  small  (39),  the  blowing 
velocity  Is  large  and  positive  as  x  ♦  xfl  and  after  peaking  becomes 
negative  In  order  to  bring  the  displacement  thickness  back  to  more  moderate 
values.  Thus  there  Is  some  parallel  here  with  the  leeslde  separation.  Further 
encouragement  to  this  notion  comes  from  the  structure,  Eq.  (39),  which  suggests 
that  In  the  central  regime  of  the  singularity,  the  velocity  components  are 
almost  constant.  If  this  Is  the  case,  In  the  leeslde  separation  singularity 
the  governing  equations  for  steady  three-dimensional  and  unsteady  two- 
dimensional  boundary-layer  flow  are  effectively  the  same  and  so  the  assumption 
of  a  van  Donmelen  and  Shen  singularity  Is  as  consistent  for  the  one  class  as 
the  other. 
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Figure  2.  Variation  of  skin-friction  components  and  displacement  thicknesses  with  p  and  T  for 
(a)  Streamwise  skin-friction  component. 


(b)  Circumferential  skin-friction  component. 


Variation  of  streamwise  skin-friction  component  U'  with  p  and  T  for 


Variation  of  skin-friction  components  ar.d  displacement  thicknesses  with  p  and  T  for 
(a)  Stream wise  skin-friction  component. 


.  (b)  Circumferential  skin-friction  component. 


Figure  5.  (b)  Streamwise  skin- 


